Exact Elliptic Compactons in Generalized Korteweg-DeVries Equations 
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We derive a general theorem relating the energy, momentum and velocity of any solitary wave solution of the 
generalized KdV equation which enables us to relate the amplitude, width, and momentum to the velocity of 
these solutions. We obtain the general condition for linear and Lyapunov stability. We then obtain a two param- 
eter family of exact solutions to these equations which include elliptic and hyper-elliptic compacton solutions. 
For this general family we explicitly verify both the theorem and the stability criteria. 



In a variety of physical contexts one finds nonlinear field 
equations for which a wide class of solitary wave solutions 
can exist. However, in many cases one is not able to obtain the 
solution in a closed form and thus it is not very easy to study 
the stability of such solutions. In this Letter we derive a gen- 
eral theorem relating the energy, momentum and velocity of 
any solitary wave of a generalized Korteweg-De Vries (KdV) 
equation of Cooper, Shephard and Sodano (CSS)i. The im- 
portant point is that this particular generalization of the KdV 
equation is derivable from an Action Principle. Using the the- 
orem, we are able to relate the amplitude, width and velocity 
of any of the solitary wave solution even if such a solution is 
not known in a closed form and also study its stability. Sec- 
ondly, we obtain a two-parameter family of solutions to these 
equations which include elliptic and hyper-elliptic compacton 
solutions. For this general family we explicitly verify the the- 
orem as well as the stability criteria. 

Compactons were discovered originally in an extension of 
the KdV equation by Rosenau and Hyman in Ref. Com- 
pactons are fundamental excitations (i.e., solitary waves with 
a compact support) of KdV-like equations that possess a non- 
linear dispersion and collide quasi-elastically. They play an 
important role in pattern formation and emergence of nonlin- 
ear structures in physical systems. Other physical contexts 
in which compactons are relevant include fluid dynamics, op- 
tical waveguides and the field of intrinsic localized modesi 
Breather compacton^ and compactons in other nonlinear 
dispersive equations^ are also known. In addition, KdV equa- 
tions with higher power of nonlinearity and dispersion may 
lead to the phenomema of blow-up and collapse. Therefore, 
the study of compactons and their stability in this class of 
equations is important in its own right. 

Rosenau and Hyman showed that in a particular generaliza- 
tion of the KdV equation, defined by parameters (m, n) (with 
m, n integers), namely 



0. 



(1) 



a new form of solitary wave with compact support is a solution 
of this equation. For the case m — n (m integer) these com- 
pactons had the property that the width was independent of the 
amplitude. In Ref. [jit was stated that K{3, 2) had an elliptic 
function solution. In a later work, Rosenai;^ obtained ellip- 
tic function compactons for the case of K{A, 2) and K{5, 3). 
Phase compactons have also been investigated^ Because the 



equations of Rosenau and Hyman were not equivalent to a 
Hamiltonian dynamical system, CSSi considered instead a re- 
lated generalization of the KdV equation 



K*{l,p) : ut + UxU 
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-r a[2Uxxx'^ 

= 0. 
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Equation (|2} has the property that it is derivable from a La- 
grangian. CSS showed that the two equations have the same 
class of solitary wave solutions when I — m+1 and p = n — 1. 
Because of this connection, the two parameter family of new 
solutions that we will find here will also be solutions of the 
K{m, n) equation with slightly different coefficients. Cooper 
and KharsS. later showed that these equations with I = p + 2, 
i.e. m = n < 3 and m continuous, had compacton solutions 
of the form 



u{x,t) ^ A[cos{l3y)f/P , y ^ x - ct , 



(3) 



and for all these compactons also the width was independent 
of the amplitude. 

The CSS equation has three conserved quantities, i.e. the 
Hamiltonian H, momentum P and mass M given by 



H = 



auP{uxy 



dx . 



P 



- I u''{x,t)dx; M 



u{x, t) dx . 



(4) 



(5) 



On the other hand, the Rosenau-Hyman equation has in gen- 
eral two conserved quantities (except when m = n in which 
case there are four conservation lawsSilli) given by 



M 



u{x, t) dx; 



u''+\x,t)dx. 



(6) 



The fact that the CSS equations were derivable from an Ac- 
tion Principle allowed CSS to consider time-dependent vari- 
ational approximations based on simple time-dependent trial 
functions u^{x,t). By using post-gaussian trial wave func- 
tions u"{x,t) = A{t) cxp[—f3{t)\x + q{t)\''] it was shown 
that these trial wave functions satisfy the relationship 



H 



(7) 



2 



with 



p + l + 2 

'p + 6-l 



(8) 



Here q is the velocity of the compacton. However, it was not 
known at the time if the resuh was true in general. In this Let- 
ter we will show that this relationship is entirely general and 
does not depend on any trial wave function approximation. 
We will show explicity that this relationship can be derived for 
all solitary wave solutions of the CSS equations of the form 
u{x, t) = AZ[f3{x + q{t))]. This relationship will also allow 
us to relate the amplitude, width and velocity of the solitary 
wave. Having a Hamiltonian formulation also simplifies the 
discussion of stability, and using general arguments we shall 
prove for the CSS equations that the compacton solutions are 
stable provided 2 < ^ < p + 6. 

Energy-momentum relation theorem: To derive the rela- 
tionship between the conserved energy and the conserved mo- 
mentum, the starting point is the action 



r = Ldt 



(9) 



where L is given by (note (f)^ = u) 



(10) 

If we assume the exact solitary wave solution is of the generic 
form 



cj),=u = AZ{p{x + q{t))), 



(11) 



then the value of the Hamiltonian for the solitary wave solu- 
tion can be shown to be 



H = -Ci(0 ^^(fl^) + a(5A^+^C2{p). (12) 

where 

Ci(0= [ Z\z)dz- C2{p) = l[Z'{z)fZP{z)dz. (13) 



Since H and P are conserved, therefore we can rewrite the pa- 
rameter A in terms of the conserved momentum P and obtain 



P^\j dxu^ = ^Cs , - y dzZ^iz) . (14) 



Replacing A by P, we now have 

H = -C^{l)P"^p'^'-^^'^ + Ci{p)P^'P+^'^'^p'^P+^^'\ (15) 
where 



Ciil) 



2 



1/2 



; Ci = aC2{p) 



2 



(p+2)/2 



(16) 



The exact solutions have the property that they are the func- 
tions of {3 that minimize the Hamiltonian with respect to (3. 
Explicit examples showing this for both the CSS equation as 
well as a quintic generalization of this equation are found in 
the appendix of Ref.[nl On using dH /d(3 ~ 0, we obtain 



/3 = PT 



C4 P + 4 



2/(/-p-6) 



This leads to 



li = f[l,p)p- 
where r is given by Eq. Q, and 



p-l + 2 
p + 4 



C4{p)p + 4 



(17) 



(18) 



(i-2)/(;-p-6) 



(19) 



Hamilton's equation, q = dH/ dP now yields the relationship 
as given by Eqs. and (8). 

Using Eqs. Q, MAI and Mil to M% it is easy to show that 
the momentum P, amplitude A and the width parameter (3 
functionally depend on the velocity c (note c = —(7) by 



P + 6-i _J 

P (X C^Ti^ , ^ (X c'- 



Several comments are in order. 



(3 (X c^rfey . (20) 



1 . Notice that when I ~ p + 2 then (3 is independent of the 
velocity c and momentum P and hence the amplitude 
A of the solitary wave. 

2. Note that the c dependence of the amplitude A solely 
depends on the parameter I and is independent of the 
parameter p. 

3. In the special case when p = and I = d + 2, the 
CSS equation reduces to d-th order KdV equation. In 
particular, in that case d — 1 corresponds to the KdV 
equation while d ~ 2 corresponds to the modified 
Korteweg-De Vries (mKdV) equation. For this case, 
a well known exact solution is f{y) = Ascch^^'^{(3y), 
where f3 = {d/2)^/v and 2A'^ = (d + 1) (d + 2)u which 
indeed is consistent with relation (20). For that case we 
notice from above that for any d, the width parameter (3 
varies as c^/^ while the amphtude A varies as c^/'^. 

Stability of Solutions. The stability problem at / = p + 
2 was studied in Ref. IT3 using the results of Karpmani^ii^. 
Their analysis is in fact also valid for arbitrary I, p. The result 
of detailed analysis is that the criterion for linear stability is 
equivalent to the condition. 



dP 
dc 



> 0. 



(21) 



Since for all of our solutions P cx c'^P+^-'^y^^^-^^ it imme- 
diately follows that the solutions are stable provide d 2 < I < 
p+6. Analysis of Lyapunov stability following Refs. Il0ll2ll3l 
also leads to the same restriction on p. 
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Exact solitary wave and compacton solutions: If we as- 
sume a solution to ^ in the form of a travelling wave: 

u{x,t) = f{y)^ f{x~ct), (22) 

we then obtain 
Integrating twice we obtain: 



(23) 



1(1-1) 



af^fP = Xi/ + K2 



(24) 



For compactons, the integration constants, Ki and K2 are 
zero. The general theorem derived above is valid in the 
case that the integration constants Ki , K2 are zero. Un- 
less stated otherwise, throughout this paper we shall consider 
the case when Ki , K2 are both zero. On demanding that 
f'F 0, f^p-^ -> at edges where / ^ 0, while 
/' is finite at edges gives us the following bounds on I and 
p : / > 1,0 < p < 2,p < L 

It is worth noting that Eq. (I24> is very similar to the equa- 
tion obtained for the Rosenau-Hyman case: 



(/') 



l\2 



2v 



n(n + 1) 



3~n 



n{n + m) 



m—n+2 



(25) 



Thus, we see that with I = rn+\ and p = n — 1 the equations 
for finding solutions are identical in form, with only differing 
coefficients. Therefore we expect to find similar solutions to 
the two sets of equations. 

Let us now look at the various different compacton solu- 
tions to Eq. MA\ . For the particular case I = p + 2 (m ~ n). 
Cooper and Khara^ were able to show that the CSS equation 
has solutions of the form (|3} and that for all these compacton 
solutions the width is independent of the amplitude. We now 
show that if instead I = 2p + 2 then one gets a one parame- 
ter family of elliptic compacton solutions. In particular, it is 
easily shown that 



for 



/ = Acn^(/3y,A:2 = l/2), (26) 



-K{k^ = l/2)< Py < K{k^ = 1/2) , (27) 



and zero elsewhere is an exact elliptic compacton solutions to 
the field Eq. (23) provided 



7 = 



2/p, l^2p + 2, A^P ^ c{p + l){2p I) , 
4 

cp 



16a2(p+ l)(2p+ 1) 



(28) 



The (m, n) — (3, 2) and (m, n) = (5, 3) are two special cases 
of solution ( I26> with 7 = 1,2. Here cn(y, k) is a Jacobi ellip- 
tic function and K{k) denotes the complete elliptic integral of 
the first kind with modulus k. 

If instead I — ip+2, then we obtain a one parameter family 
of elliptic solutions of the form 



f = A 



1 -cii(2(3)i/4/32;,fc2) 



(1 + V3) + (\/3 - l)cn(2(3)i/4/3y, fc2) 



(29) 



with the modulus 



1 

^ ~ 2 4 



(30) 



Here 



7 



1/p, ; = 3p+2, 2£'P = c{ip+\){?,p + 2), 



c^p^ 



256a3(3p+ l)(3p + 2) ' 
This is a compacton solution in the range 



(31) 



Q<w = 2(3)1/4^ < 4if |^fc2 = i _ ^ j . (32) 

and zero elsewhere, where ^ = (3y. The way this solution is 
obtained is by starting with the ansatz / = AZ°'{S) and de- 
manding that Z satisfies the differential equation [dZ /d^)"^ = 
1 — Z^. The above integral can be evaluated by converting the 
differential equation into the standard differential equation for 
the Weierstrass elliptic functionii 7^(1/, /c). On simplifying, 
we obtain the explicit solution as given by Eq. (28). The 
(jTi, n) — (4, 2) solution of Ref. □ is a special case of our 
general solution ( I29l l with p ~ 1. 

We now show that all the above solutions are in fact special 
cases of the two parameter (a, t) family of solutions obtained 
by assuming 



(33) 



and demanding that 



{Z'f = l-Z^ 



(34) 



We find the conditions: 

l=pt + 2: a^2/p, 2AP* = {pt + 2){pt + l)c , 



0- 



2t 



c'-^p' 



*23*-i(pi + l)(pt + 2) 



(35) 



It may be noted that here l,p,t are all continuous (i.e. real) 
parameters. The various solutions discussed above correspond 
to i = 1,2,3. In fact an explicit compacton solution can also 
be found in the case t — 3/2. The solution for t = 3/2 is 
essentially the same as that for t = 3, Eq. (29), but with 
Eq. (30) replaced by the complementary modulus k'^ — ^ + 

For t > 3, the compacton solutions are related to hyper- 
elliptic functionsii 

It is interesting to note that even though none of the hyper- 
elliptic solutions can be obtained in an analytic form (t > 3), 
still their momentum and energy can be obtained analytically. 
In particular, on using Eq. in the expressions for H, P as 
given by Eqs. (4) and (|5jl it is easily shown that 



H 



A^c{<6+p-l) 
'2pt {l+p + 2) 



B 



2pt 



4 1 
2 



2pt' 



B 



4 1 



2pt ' 2 



(36) 



(37) 
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where B denotes the Beta functionii. For all these solu- 
tions with p, I continuous (i.e. real) variables, the relationship 
H/P = c/r is always satisfied. 

For the hyper-elliptic compacton solutions we have that I, 
p and t are related by / = pt + 2 so that the general stability 
criterion can be rewritten as < < p + 4, or for t > 1 



0<p{t-l) < 4. 



(38) 



The requirement for non-singular solutions is that < p < 2. 
This means for t < 3, compactons with arbitrary p in the 
allowed range are linearly stable while when t > 3 , the com- 
pactons are stable only for < p < A/jt — 1). Analysis of 
Lyapunov stability following Refs. Il0ll2ll3l also leads to the 
same restrictions on p. 

From the stability analysis we also recover the well known 
result that the higher order KdV equations, characterized by 
p = 0,1 = d + 2, have stable soli ton solutions provided d < 4 
while one has unstable soliton solutions in the casai^ti^iii d > 
4. 

General remarks. In conclusion, we have obtained expUcit 
exact compacton solutions in terms of Jacobi elliptic functions 
which exist at particular values of the elliptic modulus k. Note 
that at fc = these solutions do not exist. In addition, we de- 
rived a quite general theorem that relates the energy momen- 
tum and velocity of solitary wave solutions without knowing 
the explicit form of the solution. We note that the above analy- 
sis should also hold for the generalized quintic KdV equation. 
We also notice that the radial part of the generalized nonlinear 
Schrodinger equation with nonlinear power k 



(39) 



obeys the same equation as (24) with p = and / = 2k + 2, 
so that the stability analysis of these two problems are related. 
We believe that it should be possible to derive a similar theo- 
rem in other nonlinear systems. 

Finally, as a byproduct of our results, we are able to obtain 
analytic expressions for /v (fc) and the complete elliptic inte- 
gral of the second kind, E{k), at fc^ = 1/2 — \/3/4 (see Eq. 
30) or at fc = sin(7r/12). Specifically, 



3i/40Fr(l/6) 
6r(2/3) 



(40) 



where T{a) denotes the Gamma functior^ii. In addition, at 
fc = sin(7r/12) using the relations K' = and 



E = 



rV3 
12^ 



k'K, E' 



7lV3 



k'K' 



(41) 



we also have the explicit analytic expressions for E and the 
two complete elliptic integrals with complementary modulus 

fc' = Vl - fc2 = cos(7r/12), namely K' and E' . 
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